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Abstract 

Generalizing the well-known relations on characteristic functions on a 
plane to the case of a one-dimensional regular surface (curve) with com- 
pact support, we establish implicit equations for these functions. Intro- 
ducing an approximation, we solve the combinatorial problems and reduce 
these equations to a set of linear equations for a finite number of unknown 
functions. Imposing natural conditions, we obtain a closed system of lin- 
ear equations which can be solved for a given surface. Its solutions can 
be used to calculate the distribution of hitting probabilities for a regular 
surface with compact support. 

In order to verify the accuracy of the approximate distribution of hit- 
ting probabilities, numerical analysis is being made for a chosen surface. 



Introduction 

Search for the distribution of hitting probabilities is an old and a well-known 
problem. Consider random walk on d-dimensional lattice (in continuous case 
consider Brownian motion). Then fix a surface of interest S. Suppose that any 
random walk starts from a given point z which does not lie on S. The problem is 
to calculate the distribution Pz (x) of probabilities of first contact with points x 
of the surface S. In other words, we are looking for the probability that random 
walks from z ^ S to x G S do not touch other points y E S\{x}. Of course, the 
distribution Pz{x) depends on z and S. 

This problem has been solved exactly for some particular surfaces. For 
instance, the case of a planar surface in 2D (an ordinary straight line) is de- 
scribed in any book on probability theory (see j^], Q). Its generalization for 
d-dimensional hyperplane is also simple (for example, see the end of Section 2). 
Note that exact solutions have been found only for some particular surfaces but 
not in the general case. In the general case, the asymptotic behavior is widely 
studied, |]. 

Problems of the hitting probabilities do not only have a purely mathematical 
interest. They are important for a wide class of physical problems, in partic- 
ular, for the problems of Laplacian transfer across an interface, for instance, 

^ E-mail: Dcnis.Grcbcnkov@polytechnique.fr 
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diffusion through a membrane, electrod problems, heterogeneous catalysis, etc. 
(see g, g, @ for details). Indeed, if we are interested in diffusion through 
a semi-permeable membrane (points of this membrane can absorb or reflect 
touching particle with certain probabilities), we can write the total probability 
of absorption by a chosen point of the membrane as a sum of probabilities to 
be absorbed after 0, 1, 2, etc. reflections (rigorous formalism is described in 
(H). Here we face the task to calculate the distribution of hitting probabilities. 
Note that using this distribution solely for a planar membrane, we have recently 
obtained some important results about general characteristics of the Laplacian 
transfer across an interface, (t). To solve these problems one needs to know 
the distribution of hitting probabilities for a general surface. Here we propose 
a method to approximate the distribution of hitting probabilities for a rather 
general case in 2D. 

In the first section we introduce definitions and conditions which are required 
in what follows. In the second section we briefly describe a well-known case of 
the hitting probabilities on a horizontal axis. Main results are contained in the 
third section. Section 4 is devoted to some numerical results. In the last section 
we make conclusions and discuss possible generalizations. 



1 Definitions 

Consider a square lattice on a plane. Let us define a regular surfac^ with 
compact support S = { {x, S{x)) } by a function S{x) with integer x subject to 
the following conditions: 

1. Bijection: The function S{x) is a bijcction between the set of integer 
numbers (absiccac x) and the set of surface points ; 

2. Regularity: For any x, \S{x + 1) — S{x)\ < 1 ; 

3. Compactness: 3M : S{x) — for \x\ > M , i.e. the non-plane part of the 
surface has a finite size. In other words, function S(x) has a compact support. 
Moreover, we suppose that the surface is centered : S{zk(M — 1)) ^ 0. 

Let us briefly discuss this definition. The second condition allows to simplify 
all calculations and formulae, but it does not seem to be essential (see Section 
5). Note that this assumption can be viewed as a regularity condition for the 
surface in continuous case : S'{x) < 1. 

On the contrary, the third condition is important. It tells us that the surface 
in question is a finite "perturbation" of a planar surface (line). In other words, 
this surface is composed of two parts: a complex but compact part in the center 
with two plane "tails". Moreover, it is important that both tails lie on the same 
height (which is chosen as 0). This feature will allow to obtain an approximate 
distribution of hitting probabilities by using the same ideas as for a planar 
surface (see Section 2). 

We call all the points {{x,y) : y = n} the n*'* level. Denote 
N = max{5(a;)}, N* = - min{S'(a;)}, 

i.e. the surface lies between (— iV*)th and A^th levels. 

^ Even for two-dimensional case we prefer to use the word "surface" instead of "curve" or 
something else. 
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All points M = {{x^y) : Vx y < S{x)} arc called internal. All points 
£ = {{x,y) : V.T y > S{x)} are called external. The external points near the 
surface, {{x,S{x) + 1)} are called near-boundary points. The functions defined 
on these points, are called near-boundary functions (see below). Often we shall 
use the words "surface" , "near-boundary functions" , etc. thinking only about 
the non-trivial part, i.e. for \x\ < M. 

The external points with y ~ arc called ground points. The functions 
defined on these points, are called ground functions. Let J = {k : (/c,0) e £} 
the set of abscissae of ground points. Let also Jo = {k G (— M, M) : {k, 0) e 5} 
the set of abscissae of boundary points on zeroth level (only non-plane part!). 

We introduce the hitting probabilities Pk,n{x), i.e. the probability of the 
first contact with the surface at point {x,S{x)) if started from {k,n). Their 
characteristic functions 4>m,n{()) are 

oc 
x= — oo 

The inverse Fourier transform allows to obtain Pk^n{^)-, 

TT 

/rlG 
^e-"^0fe,„(0). (2) 



Planar surface 



At the beginning, we consider the trivial and well-known case of a planar surface 
(horizontal axis): S{x) = 0. This case is useful to remind the technique of 
manipulation with characteristic functions. 

Suppose that n > 0. The probability Pk,n{x) satisfies a simple identity 
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Pk-i-i,n{x) + -Pfe-i,n(a;) + Pfc,„+i(a;) + Pk,n-i{x) 



which can be also written for characteristic functions, 



4>k+l,n{6) + <i>k-l,n{6) + <i>k,n-^l{6) -\- (t>k,n-l{0) 



Translational invariance along the horizontal axis gives 

0fe,„(e) = e«'=Vo.n(^^). 
Using the obvious condition Pk,oix) = Sk,x, we obtain 



■Ae) 



ike 



(3) 



(4) 



(5) 



(6) 



The last trick is the following. If the starting point is placed in the n-th level, 

the random walk must cross the {n — l)-th level at some point (m, n — 1) to 
reach zeroth level. The probability to pass from {k,n) to {m,n — 1) without 
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touching other points in the (rt — l)-th level is exactly Pfc^i(m). Therefore we 
can write 

Pk,n{x) = Pfc^i (m)P„i,„-i (x) . 
m 

In terms of characteristic functions this convolution is just a product of the two 
corresponding characteristic functions, 

Using the translational invariance (^), we obtain 

^fc,„W = e'^''[0o,iW]". (7) 

Substitution of expressions (^) and into relation (^) for n = 1 and A: = 
leads to 



MO) - J ( e-'Vo.iW + e'Vo.i(^) + 1 + [00,1 W]' 



<?!)g.i-(4-2cos6l)</.o,i + l = 0. 



or 



(8) 



This quadratic equation has two solutions, and we should choose the one for 
which < 1 (property of characteristic function). It is denoted ifi{9), 



<fi{9) = 2 - cos 6* - V(2- cos6')2 - 1. 
So, we obtain for the planar surface 



(9) 



(10) 



Inverting this relation with the help of (||), we obtain the distribution of hitting 
probabilities for the planar surface, 



ryplanar 
k.n 



— I 

2tt 



(11) 



This well-known result will be used for a general case (some properties of co- 
efficients i?^ are described in Appendix 1). The formulae (H) and (|l^) can be 
generalized for d-dimentional hypercubic lattice. 



^p{9i...9d-i) =d-^ cos( 



i=l 



K,n((^i---fd-i) = exp 



d-i 



i=l 



d-^cos(0,) -1, 

i=l / 

¥'"(^i...ed-i). 
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3 Regular surface with compact support 



We shall consider the characteristic functions 4>k,n{0) as a vector 

$(")(6/) = 



(0) 



of (2L + 1) components where parameter L is supposed large, and it will tend 
to infinity at the end of calculation. 

For the planar surface we had relation (^) which can be written in matrix 
form 

= $("-!) + $("+1)^ (12) 



where the matrix A is 

Ai^i = 4, Ais+i 



A. 



I7 A2L+1,1 — Ai^2L + l 



-1. 



The last equalities are artificial : we added them to obtain a cyclic structure of 
A. But at the limit L 00 this little modification vanishes. The eigenvalues of 
A are 

Ah = 4 - 2 cos(6ift), with e,, = 



2L + 1' 



and the eigenvectors are given as 



Vh 



' g-*(L-l)0h 



V 



I 



Now we generalize the relation to the case of a regular surface with compact 
support by introducing vector A<i>*^"\ 

^$(") = $("-!) + $(»+!) + A<I>*^"^ (13) 

(this relation can be regarded as the definition of A<i>'"'). Let us introduce 

c„(^?,(?,) = ($("), K), Ac„((?,0;,) = (A$("),y^*). (14) 

We can rewrite ( p^ in terms of c„ and Ac„, 

AftCn = C„„i + C„+i + Ac„. (15) 

If we can express c„ in terms of A/j, ^{Q\ cq and {Acfc}, we find <i>*^"'' as a 
decomposition in the eigenbasis 



(16) 



(factor (2i + 1)~^ is due to normalization (V/i,VJj*)). In order to solve the 
recurrence relations (|l5|), we should close them by certain conditions. For the 
lower half plane we take a sufficiently large number TV; > iV* , and 



^-N, = 0, 



c-jv, = 0, 



(17) 
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because it is impossible to penetrate through the surface. More generally, ac- 
cording to the definition of hitting probabilities, we should maintain 

[ 0, if (to, n) G M. 

For the upper half plane we shall use the following trick. We take a large 
number ^ N and consider 4>x,Nu+i{^)- fo'" ^ planar surface, we can 
write Pa;,Ar„+i(ri) in convolution form. 



where H]^_^ is the probability to hit point (to, N^) if started from (a;, 
For characteristic functions it is simply 



In the plane case we used the translational invariance along the horizontal axis 
to simplify this sum. Evidently, such a symmetry breaks down in the general 
case. But if we take Nu sufficiently large, i.e. we "look" on the membrane from 
a remoted point, we can suppose that translational invariance is approximately 
valid, 

'/>m,iv„(^)«e'^™"")Vx,iv„(^). (19) 
Using this approximation, we immediately obtain 



(0) ^ ipie)^xMA0), 



where function ip is defined by (|^) . This is our approximation in the upper half 
plane which allows to close the recurrence relations (|l5|), 

CN^+i{e,eh) = ipie)cNje,eh). (20) 

The main idea is to step down from A'^u-th and (— A^i)-th levels to zeroth 
level. We shall consider the upper and lower half planes separately because 
the relations (|l^) and ( pO| ) are different. Note the essential complication of 
the general case with respect to the planar surface. For the planar surface 
we had Ac„ = for any n, and the system of equations ( |l5| ) was closed. It 
was sufhcient to solve these recurrence relations by substitution c„ — cqc"" in 
(p^, and we obtained the final form of On the contrary, for the general 

surface Ac„ ^ 0, and they depend on the near-boundary functions (j)m,n (see 
below). Consequently, the decomposition ( p^ ) itself becomes a system of implicit 
equations for 4>m.n- For the moment, the problem is complex. It will be solved 
by two steps. First, we obviate the combinatorial problems, i.e. we find the 
explicit solution of reccurence relations (|l5|). Second, we solve the equations for 



3.1 Solution of reccurence relations. 

A direct verification shows that 

Cn = Pn^-iiCN^ - ^ a;Ac„+i, (21) 
1=1 
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Afj-n-l 

C-n = aAr,_„C_jv, + l - ^ a;Ac_„_i (22) 

is a general solution of (|l^) (we omitted the index h which does not change the 
structure of the solution), where 

ao = 0, ai = 1, an+2{0h) = '>^han+i{Oh) - ctniOh), (23) 

Pn{e,eH) = an{eh)[Xh - ^{e)] - a^-i{9h), Po = i- (24) 

We are looking for the explicit representation for q;„ in the form 
Substituting this into (|2|), we obtain equation 

— XhX + 1 = 0, 

which has two well-known solutions: ip and (compare this equation with 
(H)). As the expression (|2^) is linear, we find a general solution as linear com- 
bination of and such that ao — 0. We obtain 

««(M = \^^j^</.^-"(e.), (25) 
or as a geometrical sequence 

n-l 

a„(^,) = ^¥>2^'+i-"(^/0- (26) 
i=o 

Formulae ( pT| ) and ( p^ are valid for any n > 0, in particular, for n — and 
we can express cm^ and c_jv,+i in terms of cq and {Ac;}, 

1 / ^" \ 1 / \ 

= ^5 — Co + a; Aq , c_Ar,+i = + a; Ac_i , 



hence 

co + ^aiAcij- ^ a/Ac„+i, (27) 



OI-Ni-1 



1=1 / 1=1 

Ni-1 \ Ar,-n-l 

C0+ I - ^ a/Ac-„-i. (28) 



' \ 1=1 / 1=1 



Let us introduce 



PN^ OLNi 

(the right-hand side of fi^"^ depends on 9 through factor ip; the dependence on 
9fi is due to a„ which contains A/i). 

Using the reccurence properties of a„ (see Appendix 2), we simplify relations 
(13) and (|2|), 

/ n \ iV„ 

c„ = /(^"Mco-l-^aiAQ -!-«„ ^ //^"^Aq, (29) 

\ i=l / l=n+l 

/ n \ W,-! 

c^n = /r'Mco+^a,Ac_, +a„ ^ //^')Ac_,. (30) 

\ 1 = 1 I l=nJrl 
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Rewriting the definition (|2J) of /?„ as 
after simplifications we obtain 

We remind that is an arbitrary sufficiently large number. Therefore we can 
take the limit oo. Knowing that (p{0h) < 1 for 9h ^ 0, we obtain that in 

this limit (^2Ar„-n+i ^^^^ ^2N^+i vanish, i.e. 

/(-)(0,(?,O = ^"(e/.)- (31) 
Using formula (|25|), we can also write the explicit representation for fn^'^\ 

As above, we take the limit A^j — > oo to obtain 

fi^\0H) = ^'\9h). (32) 

Note that we obtained the same limits and imposing the different 

conditions (|l^ ) and (po|) in the lower and upper half planes respectively. Using 
formulae (|3l[) and (|32|), we can write 

Cn \ oc 

co(0,0/,)+^ai(^/OAQ(0,0/.) +«„(0/,.) ^\Sh)Aciie,eh), 
1=1 ) 

(n \ cxD 

Co(^,0/.)+5I"K^/OAc_,(f?,0,.) +a„(0,.) ^'iOh)Ac_i{0,0h) 

Introducing functions 

( Lp'^{e)ai{e), if / < n, 

l'l''\e)^\ ip'{0)a^{0), a l>n, (33) 

[ 0, if / < or 77. < 

(the last convention will be used in the following), we can write c„ in the unique 
form (for 77 > and 77 < 0), 

N 

c„ = ^i"i(0;,)co(0,0,)+ Y bl''\o,,) + j[-''\0^)]Aci{e,e,,), (34) 

where summation over I is just from —N* to N, because corrections Acn+i and 
Ac_jv*-i are equal to for /> (see Section 3.4 for details). 
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3.2 CoefRcient cq 

Let us calculate the coefficient cq, 



fc=-L 



As it was mentioned above, the plane "tails" of the surface lie on the zeroth 
level, thus 

(<i>o)fe = e^'=^ if |fc|>Af. 

Using this explicit form, we are going to compute the contribution cj^ ^ of plane 
"tails" 

(0),. „^ ^ , _ sm{L + O.5)i9-0h) sin(Af - 0.5)(g - g,,) 



4°^(0,0,O = 2 E cos(0-0,O = 



^^^j sin(0 - 9h)/2 sin(0 - 0ft.)/2 

(35) 

(in the case 6 = 9h one should consider this relation in the limit sense when 
^ 9h)- Note the simple relation 

dO .xSin(2M-l)0 r 1, if kGA, 

-cos{2ke) — =X(-M,M)(fc), = if^^^^ 

(36) 

The contribution Cq^-* of intermediate part of the surface (with \x\ < M) can 
contain some non-trivial functions (t>m,o{9), with /c G J and k £ Jq 

ci'\0,eh) = 5] + E e^'^'"''^'- (37) 



3.3 Limit L — > cxd 

Expression (^) transforms to the integral expression for (/)/j_„ by taking the limit 
L oo. Here we write only the expression for n > 0, the opposite case will be 
easily obtained later. Note that the first term of (|35| ) tends to ^-function in the 
limit L — > oo. It removes the integration in the first term, i.e. we obtain 

'^fe,«(^) = ^"(^)e^''-C(^)+ (38) 

^e^'"' (^^"(^')c^'^(^?,^')+E7,'"\^^')Aq(0,0')) • 

The first term is the contribution of plane "tails" . The second term corresponds 
to the perturbation on zeroth level due to Cg*^^ , 



The last term of ( pSf ) contains some unknown functions (pm.iiS) through the 
coefficients Aci{9, 9'). It is denoted as T[0], and we are going to calculate it. 
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3.4 CoefRcients Aq 

To get ahead with the expression (^ ) , we should write expHcitly the coefficients 
Aci{9, 9'). It is not so easy for the general case. Indeed, for these purposes one 
calculates the contributions of each point on l-ih level. The problem is that there 
are many conditions, and they lead to complex formulae difhcult to manipulate. 
We are going to present the other way. 

What is the origin of the vector A$(")? Let us recall the definition (hSf) 
where vectors A^^"-* were introduced to generalize the expression (|l|). A brief 
reflection shows that 



(A$ 



0, if (m, n) e 



In other words, the relation ( |l^ ) is satisfied automatically for any external point, 
but it should be imposed artificially for each surface and internal points. 

Now it is the moment to remind the formula (|l^ which tells that functions 
4'm,n{9) are equal to zero on the internal points. Therefore we can consider only 
the points near the surface S. A direct verification shows that 

r e^(™+i)«, if S{m + 1) - S'(m) = 

(A<i>(^(™)))„ = 46™*^ - 0™,5(™) + l - I <t>m+l,S(rn), H S{m + 1) - S (m) = -I 

[ 0, if S{m+l) - S{m) = 1 

g»(m~i)e^ if S{m - 1) - S{m) = 
(f>7n-i,s(7n), if S{m - 1) - S{m) = -1 , (40) 
0, if 5(m- 1) - S{m) = 1 

(A?>^ 'jm — -e -e' Os{ra),S{m+l) + l - e <^ S(m) ,S{m-l) + l, 

(A$(")),^ = 0, if n ^ S{m) and n ^ S{m) - 1 

(here we use Kronecker 5-symbols, Su — 1 and dij = if i 7^ j). Usually there 
are several nonzero components of A^^^"^ for each n, because each level contains 
several surface points. But there exists one exception - zeroth level, where there 
is infinity of surface points due to the plane "tails". Thus, the vector A^f^"^^ 
has exceptional structure^. It contains the usual terms due to the non-trivial 
part of the surface, and the contribution of plane "tails" . Note that the last one 
is equal to —c^"^ which was calculated in Section 3.2. Later we shall use this 
result for the lower half plane. 



3.5 Approximate distribution of hitting probabilities 

According to the definition (|^ of Aq, we can rewrite T[0] as 

_^ m=-M+l 1=1 

^ There exists another vector with exceptional structure, A$(''). However, jap construct 
our treatment by the way when this vector does not appear in expressions (see {B9|), (pd)). In 
other words, we step down from A^^-th and (— Af;)-th levels to zeroth level in order to do not 
pass through the zeroth level. 
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where we changed the order of summation over m and I. However, in the last 
sum there are only two terms corresponding to (A$'^'^(™)^)m and (A<i>*^'^('"^~^^)m, 
if S{m) > 1 (otherwise, this sum is equal to 0). Using expression we obtain 
explicitly 



tn=-A/+l. 



2^' 



4e -<Pm.,S(m) + l-e OS(m),S{m+l) 



''l^m+l,S(m)5 S(m),S(m+l) + l ~ 6*^™ "^^^<5s(m),S(m-l) ~ 0m-l,S(m)<5s(m),S(m-l) + l 



+ ' Os(m),S(m~l) + l+ <i S{m) ,S(m+l) + l 



(here we have used the last convention in the definition (^) of 7}"' to avoid 
any terms with S{m) < 0). 

The last step is to transform this huge expression for characteristic functions 
into hitting probabilities using the formula (||). Note that all functions e*™^ after 
integration over 9 with e"*'^^ give J-symbols that remove the summation over 
m in corresponding terms (but we should write factor X(-Af,Af) (2^))j 



T[P] 



M-l 
m=-A/+l 



(41) 



(m).S(m+l) + l ~^ ^rn— l,5(m) (m).S(m— 1) + 1 I f 

Xi-M,M)\^) [^^x,S{x) - ^x-l,S(x-lfS{x),S{x-l} - ^^+i^s{x+l)^S(x),S{x+l)- 

■p.(k,n) T^{k,n) r ri'^k,n) ^ \ 

^x,S{x)-l - iJx+l,S{x+l)~l^S(x),S{x+l)~l - i^:r-l,S(s^-l)-l'^S(a;).S(a;-l)-l I, 



where 



(42) 



Note that coefficients d'"^'^^ are universal, they do not depend on a given surface. 
It means that once calculated, these coefficients can be used for any hitting 
problem in 2D. They can be also expressed in terms of H^, 



D 



min{n,/} 



2j-l + \l-n\ 



(if n or Z is equal to 0, the sum is also equal to 0). We just indicate several 
useful properties of these coefficients. 



p.(A:,n) pj(0,n) j-^(k—rn.n) j-^(m—k,n) 

^m,l — ^m-k,l — ^0,1 ~ ^0,1 ' 



j^(k,0) _ j^{k,n) _ „ 



j~,{k,n) j~,{k,n) r^ik.n) j~,{k,n) _ . j~,{k,n) r r 
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The first part of (^ij), containing Pm.s(m)+i, can be represented as 

M 



m=-M 

with coefficients 



The second part of (41) can be simplified. Indeed, using the properties of 
5-symbols, we have 



To = 



Xi-M,M)i^) [^^x,S{x) ~ ^2;-l,S(x) " ^x+l,S{x) ~ ^x,S(x)-l 
+^x-l,Six)'^Six).Six-l) + l + ^x+l,S(x)^Six),S{x+l) + l I ■ 



Using the properties of dI^'^\ we finally obtain 



/ \ / j-^{k,n) „(fc,n) r „(fc,n) r \ 

X(-M,M) (a;) I i^^,s(x) + l + -'^x-l,S(x)'^.SM,S(i.-l) + l+i^^+l,5(^)'5s(a:),S(x+l) + l I • 

(43) 

Here we omitted term Sk,xSn,s(x) supposing that starting point (fc, n) does not 
lie on the surface. 

Let us get back to the formula ( |3^ ) . Using the inverse Fourier transform (H) , 
we write 

Pfe,„(x) = - Pt°l{x) + J2 H'k-m5m,x + Hk_„,PmA^) + Tl[P] + T^. 

(44) 

The second term is 

p(o). ^ M^' „,^„sin(A/-O.5)(0-0') 



27r y 27r ' ' sin(6'-6")/2 

— TT — TT 

Replacing in the first integral 9i = 6 — 9', we factorize these integrals. The 
first factor is exactly IIJ!_^. The second one was calculated explicitly, see (|36|), 



and it equals to XI-m+i,m~i]{x). Consequently, the first two terms in (|44|) can 
be grouped into HJ!_^X{-oo,-M]n[M,+oc){x)- It means that the solution H^_^ 
of the plane case is valid only for the plane "tails" , whereas on the non-trivial 
surface (for \x\ < M) the main contribution is due to other terms. So, we have 
obtained an important result, 

M 



PkA^) = Pk,nix) + J2 HJ!_„^P„M - G(,^"^P™,s(™)+i(a:), (45) 

ra^J m— — M 

3re 
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(the third term is due to the second sum in (37)). Using (p3|), we can combine 
first two terms to obtain 

for 71 > 0: 

Fk,n\X) - i>'^.5.(^) + i+i>'^_i.5(^)ds(x),S(x-l) + l + i^^+l,5(^)ds(a;),5(a:+l) + l+-«fc-xXJ^ 

(46) 

In order to obtain analogous results for the case n < 0, we remind that 
expression ( p^ contains two terms: 7^'"^ and 7i;"''- In the previous treatment 
we have used only the first term. It means that analogous results for n < can 
be easily obtained by "reflection" of all "ordinates" with respect to horizontal 
axis, 

for n < 0: 

M 

PkA^)^PlM)+T.Hk-mPmA^)- E G(^")*P„,5M+i(x), (47) 

m£j m— — Al 

with 

Kni^) = X(-MM) i^) [d^x.-sU-I + ^i-W(.)^S(x),S(x-l)-l+ (48) 
'^ni - -S(m)+''S(m),S(m+l) + l-t^,„+i _s(,„+i)+''s(m),S(m-l) + l-t^„j_l _S(m^ 

In (^) there appears function X{-m,m){x)j because for n < there is no con- 
tribution Pjfnix) due to plane "tails" (see remark at the end of Section 3.4). 

Note that we cannot represent analogous expressions ( ^ ) and (^Sj) uniquely 
by writing |n| and |S'(a;)|. It is due to the fact that functions Pm,n hi the upper 
half plane (n > 0) have no influence on functions Pm.n in the lower half plane 
(n < 0) (except through the ground functions), and vice versa. For example, in 
the sum of near-boundary functions (the last term in (H) and (|3)) coefficients 
Gm'"'' should be equal to if n > and S{x) < or if n < and S{x) > 0. 

Pk.n{x) can be considered as first approximation to Pk^n{x). Note that a 
priori there is no reason to neglect the second and the third terms in (Ea) . 
Normally, we should take these terms into account, thus the relation ([45D is 
considered as a system of linear equations on the near-boundary and ground 
functions. In order to find these functions, we should close the system of linear 
equation. Note that equations ( ^5|) must be satisfied for any fc, n and x, and 
we can choose appropriate values of k and n. To close the system for near- 
boundary functions, we take {(fc,7i) : k e [— M, Af], n = S{k) + 1}, i.e. for any 
k e [-M, M] 

M 

m^J 771— — M 

(49) 

To close the system for ground functions, we can choose different conditions. 
For example, if we consider surfaces with S{x) > on a; € {—M, M), there are 
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no ground functions, thus there is no additional condition other than (|49|). For 
the general case (where J 7^ 0) , we propose the following condition 

PkA^) = ^(^Pk+i.o{^) + Pk-iA^) + PkA^) + Pk-i(x)^ ^ keJ. (50) 

We started from this relation for all the external points (see (^). Here we just 
demand that this relation remains valid if we substitute our approximations for 

Pk,i{x) and Pk-i{x). 



4 Some numerical verifications 

In this section we briefly present some numerical results to check the validity of 
the approximation. 

First of all, in Fig.l we depict the function (p{9) which plays a central role 
in this work. 




0.1 1 ^ ^ ^ ^ ^ ^ — I 

1 2 3 4 5 6 

e 

Figure 1. Function f{0). It changes in range from (p{0) = 1 to (p{Tr) = 3 — \/8- 

We can see that (p{6) <C 1 if 61 is not in vicinity of 27rm (m e Z). This 
property was used in the limits (^) and (|3^) . 

Let us consider again the planar surface, S{x) = 0. Without simulations, we 
easily obtain that 

ft,„(x) = PkMx) = 
So, in this trivial case our approximation gives the exact result (cf. (pj])). 

For a particular non-trivial surface the accuracy of the formula ( ^5|) can be 
obtained by comparing its values with numerical simulations of random walks. 
We have taken a simple surface represented in Fig. 2. 
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-10 10 



Figure 2. A simple surface with : N — A, N* — 0, Af = 10. 

In this case there are no ground functions. On the contrary, there are 21 near- 
boundary functions which can be calculated with the help of (^). We present 
two distributions of hitting probabilities obtained numerically and through for- 
mula @ (see Fig.3). 

We can conclude that our approximation is quite good. 





Figure 3. Distributions of hitting probabilities (in log-scale): the probability of 
the first contact with point (x, S{x)) of the surface if started from the point (0, 5) 
(Fig. 3a) or (15,1) (Fig. 3b). 
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5 Conclusions and possible generalizations 



Let us sum up what has been done. Using the same technique as for a planar 
surface, we obtain the reccurence relations for coefficients c„. For the lower 
half plane we impose c_jv, — for a sufficiently large iV;. This condition tells 
that random walks cannot penetrate through the surface. For the upper half 
plane there is no such condition. However, for a sufficiently large we can use 
an approximate condition cn^+i — y^CN^ supposing that from a remoted point 
the regular surface with compact support looks like a translationally invariant 
object. Then we find the explicit solution for reccurence relations under these 
conditions. The influence of our approximation becomes more and more neg- 
ligible with increasing of Nu- Taking the limits Nu — > oo and Ni oo, we 
express c„ in terms of explicit functions if and 7}"'' and coefficients {Ac;}. Note 
that two limits /n""' and fn°°'^ for the upper and lower half planes are identical. 
It means that in the upper half plane we could use the condition cn^+i = 
instead of cat^+i — ipcN^- In other words, one could imagine an absorbing line 
y = Nu + 1, and then send it to infinity (limit Nu — s- 00). These two possible 
approximations give the same final result (pT]). 

Having made these calculations, we express Ac„ in terms of near-boundary 
functions and combinations of exponential functions. Finally, we obtain a sys- 
tem of linear equations ( |4S| ) for near-boundary functions Pm.s{m)+i{x) and 
ground functions Pm.oi^)- K can be solved, and after that one can use the 
approximation ( p5| ) for any point (k,n). 

Numerical analysis shows that this approximation works quite good. 

The main conclusion is that we have found an approximate distribution of 
hitting probabilities for a rather general surface, pending certain conditions. In 
particular, one can make use of these results for a further study of the Laplacian 
transfer problems. 

Now one needs to study the role of conditions which were imposed in the first 
section. As we said above, the compactness condition is the most important. It 
tells that 

- the surface has a compact support, i.e. there is only finite "perturbation" 
of the planar surface ; 

- the plane "tails" have the same height (zero of the vectical axis). 

If we want to consider a surface with infinite support, we can obtain the 
same results but with an infinity of near-boundary functions. Thus the system 
( ^9| ) has infinitely many equations, and we cannot proceed any further. The 
same difficulty appears if the plane "tails" have different height: while we step 
down from the N^-th level to the level of a lower "tail" , we must pass through 
the level of a higher "tail". It means that there appears again an infinity of 
near-boundary functions. Only if we step down from the iV„-th level to the 
zeroth level (and from the (— iV/)-th level to the zeroth level), we can avoid the 
appearance of an infinity of near-boundary functions. 

The regularity condition is used to simplify certain expressions. Neverthe- 
less, it does not seem to be restrictive. Normally, to apply the technique of 
characteristic functions, we enumerate all sites (points) of the surface. To sim- 
plify the problem, one can make one of the following assumptions: 

- either suppose that the surface obeys the regularity condition; 

- or be interested in the total hitting probability Pjf^*^"'^^ (x) of points {x, S{x)), 
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{x, S{x) — 1), (x, S{x— 1) + 1) (if we authorize changes of S{x) by more than 
one unit). In other words, we could identify the points of the surface which 
have the same x-coordinate. Any of these assumptions allows to enumerate the 
points of the surface with their x coordinate using function S{x). In the first 
assumption we consider the surfaces having only one point for each x; in the 
second assumption the surfaces can have some points with the same x, but we 
are interested in the total probability for each x. 

In order to generalize the method, one can introduce another parametrization 
of the surface. One possible generalization will be presented in our forthcoming 
paper. 



6 Appendices 



6.1 Coefficients HJ} 

As we have seen, coefficients HJ! play a central role in all calculations of hitting 
probabilities in 2D. Here we briefly present some useful properties of HJ}. In 
real form (nil) becomes, 



TTll 



— coa{ke)ip''i9). 

TT 



(51) 



We can write two inequalities for 6 e (0, tt) : 



1 



cos(6i/2)' 



which can be useful for estimations. 

Now we are going to calculate the asymptotics of HJ^ for large k. Integrating 
the expression (|5j) by parts four times and using the values of the derivatives 
(f''{s) at the points and tt (see Table 1), we obtain the asymptotic behaviour 



n 



n{n? — 0.5) 



+ 0(A:-6), fc>l. 
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1 


-1 


1 


-1/2 


TT 


3-\/8 





3\/2/4 - 1 






Table 1. The values of the derivatives tp^''^{s) at the points and n. 
The formula (B2h works rather well for fc > 10. 



n\k 





1 


2 


3 


4 


5 


1 


0.3633 


0.1366 


0.0609 


0.0319 


0.0189 


0.0124 


2 


0.1803 


0.1221 


0.0756 


0.0477 


0.0315 


0.0219 


3 


0.1136 


0.0958 


0.0715 


0.0517 


0.0376 


0.0278 


4 


0.0826 


0.0759 


0.0631 


0.0501 


0.0392 


0.0307 


5 


0.0651 


0.0620 


0.0548 


0.0464 


0.0384 


0.0315 



(52) 



Table 2. The values of coefficients HJ} for small k (n in range from 1 to 5). 
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Table 2 shows values of HJ^ for small k. The asymptotics of HJ^ for large n 

is 

i.e. we obtained the same behaviour as for the brownian motion. It is quite a 
reasonable result : if we look on the surface from a remoted point, there is no 
difference between continuous and discrete cases. 

6.2 Manipulation with coefficients a„ and /3„ 

Here we present some properties of coefhcients a„ and /?„. Also we prove the 



formula (|29|). Using only the definition (23), we find 



a„ = aian^i+i - ai-iUn-i (53) 
for any I < n. Also using (p^, we have 



Let us prove (|29|). According to (27), we have 



/3. 



Co + ^ a/Aq j + - — ^ Acn+i{PN^-nan+i ~ PN^ai). 



f^^- \ i~i J z=i 
Now we should simplify the difference in brackets in the last sum. 

/3Ar„ - /3jv„a; = {aN^-n+i - '/'aw.. - {otN^+i - "^aNjai. (54) 



Consider the difference A — oiN^.-n+i'^n+i ~ ctNu+i'^i- Using the property 
we can reduce the index {n + I) in the first term and (A'„ + 1) - in the second 
term, 

A = aN^-n+i(aian+i - Q!,_iQ!„) - {an+iaN^~n+i - anaAr„-Ti)a; = 
(we used the property (|3|) in the last equality). Thus, we can represent (|^ as 

hence we find the formula (p9|). 



c 



n — Jn 



/=1 / ;=n+i 
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